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ABSTRACT 


The effects of shoaling inner shelf solitons on the bottom boundary layer have 
been observed and analyzed over a two month summer period at the Monterey Inner 
Shelf Observatory in Monterey, CA, during 2002. Utilizing CTD data to characterize the 
temperature field of the water column. Acoustic Doppler Current Profiler (ADCP) data to 
measure the velocity structure from 3m height above the bed (HAB) to the near surface, 
and Bistatic Doppler Velocity Sediment Profiler (BDVSP) data to measure the velocity 
structure and sediment concentration from a range of 60cm to 1cm HAB, solitary internal 
waves and internal tidal bores were regularly observed at the observation site. These 
events were characterized by their large isotherm displacements and the sudden change 
from near surface to near bed stratification as the internal tidal bores passed the 
observation site. Cross-shore timeseries revealed that the strongest events pushed water 
onshore near the surface and offshore near the bed, indicating a baroclinic water column 
during their passage. To analyze their effects on the bottom boundary layer, 3m HAB 
ADCP and BCDV velocities were compared with backscatter data and surface gravity 
wave energy at 3m HAB to determine their relative contribution to bed stress and 
resulting sediment suspension. As the strong internal waves passed, a logarithmic layer 
formed indicating that shear stress above the bed was occurring. This allowed the 
friction velocity within the log layer to be estimated. Combining this term with the stress 
contribution due to the wave energy, the total stress on the bed was then estimated. From 
this it was determined that typically moderate surface gravity wave forcing at the bed 
suspended sediment, while solitary internal waves and internal tidal bores continued to 
transport suspended sediment offshore near the bed. 
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Figure 6. 


Location of Monterey Inner Shelf Observatory (MISO) site. The 
Monterey Canyon can be seen bisecting the bay with its steep walls and 
complex bathymetry. The V denotes the location of MISO and the 3 
bars indicate possible generation sites for internal bores. The 2 stars 

indicate buoys Ml (inshore) and M2 (offshore). Courtesy of MBARI.3 

Wind speed and direction from MBARI Ml and M2 buoys. Data is 
averaged every 7 days over the period of 1 year. Scalar averaging makes 
winter storms much more apparent as seen by the arrows. During the 
winter period, these winter storms disrupt the normal wind patterns and 
help to lower the mixed layer in the water column. They help to mix the 
upper part of the water column and decrease the insulation. These factors 

combine to deepen the stratification during the winter.6 

(a) Temperature timeseries of Ml mooring in Monterey Bay. 13Nov01- 
17Oct02. (b) Salinity timeseries of Ml mooring in Monterey Bay. (c) 
Sigma-theta (Ge) timeseries of Ml mooring in Monterey Bay. Strong, 
shallow stratification can be seen during the summer months. The winter 
months show a weaker and deeper stratification. Conditions during the 
summer months are more ideal for the generation of an internal tide and 

solitions.8 

(a) Temperature timeseries of M2 mooring in Monterey Bay. 12Apr01- 
08May02. (b) Salinity timeseries of M2 mooring in Monterey Bay. (c) 
Sigma-theta (Ge) timeseries of M2 mooring in Monterey Bay. Strong, 
shallow stratification can be seen during the summer months. The winter 
months show a weaker and deeper stratification. Conditions during the 
summer months are more ideal for the generation of an internal tide and 

solitions.9 

Summer conditions at Ml mooring in Monterey Bay for yearday 172 to 
266. (a) Temperature profile indicates a strong and shallow stratification 
in the water column, (b) Salinity profile shows a strong and shallow 
stratification as well, (c) Sigma-theta (Ge) profile agrees well with the 
temperature and salinity profiles showing a strong and shallow 
stratification, (d) N (buoyancy frequency) calculated as a summer mean 
using 20m Ge differences gives indication that shallow stratification is 
occurring in the water column. The strong, shallow stratification in the 

summer is ideal for the generation of an internal tide and solitons.11 

Summer conditions at M2 mooring in Monterey Bay for jearday 172 to 
266. (a) Temperature profile indicates a shallow and strong stratification 
in the water column, (b) Salinity profile shows a strong and shallow 
stratification as well, (c) Sigma-theta (G@) profile agrees well with the 
temperature and salinity profiles showing a stronger and shallower 
stratification, (d) N (buoyancy frequency) calculated as a summer mean 
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using 20 m o 9 differences and plotted vs. depth confirms that the water 
column has a very strong and shallow stratification. The strong, shallow 
stratification in the summer is ideal for the generation of an internal tide 

and solitons. 12 

Figure 7. Winter conditions at Ml mooring in Monterey Bay for 2001-2002 yearday 
356 to 80. (a) Temperature profile indicates a very deep and weak 

stratification in the water column, (b) Salinity profile shows a deeper 
stratification as well, (c) Sigma-theta (ce) profile agrees well with the 
temperature and salinity profiles showing a deep and weak stratification. 

(d) N (buoyancy frequency) calculated as a winter mean using 20m o 9 
differences plotted vs. depth clearly indicates that stratification is weak 
and occurs very deep within the water column. The weak, deep 
stratification in the winter is not ideal for the generation of an internal tide 

and solitons.13 

Figure 8 . Winter conditions at M2 mooring in Monterey Bay for 2001-2002 yearday 
356 to 80. (a) Temperature profile indicates a very weak and deep 

stratification in the water column, (b) Salinity profile shows deeper 
stratification as well (c) Sigma-theta ( 09 ) profile agrees with the 
temperature and salinity profiles indicating that the water column has a 
weak and deep stratification, (c) N (buoyancy frequency) calculated as a 
winter mean using 20 m o 9 differences plotted vs. depth clearly indicates 
that stratification is present, but very deep in the column. The weak, deep 
stratification in the winter is not ideal for the generation of an internal tide 

and solitons.14 

Figure 9. CTD data taken from NAVO sponsored Calcofi cruises. Location of cast: 

36.74°N, 122.027°W. (a) July 2000. Strong, shallow stratification is 
clearly observed in the Temperature and Density Anomaly vertical 
profiles. The mixed layer can be observed at l-2m depth. N vs. depth 
shows Nmax occurring at the surface (l-2m) with a magnitude of 0.03 s' 1 . 

(b) August 2001. During late summer, stratification starts to move deeper 
into the water column. It can be found at 10m depth. Nmax is also seen 
to decrease to 0.015 s ' 1 and move down as well (10m). Summer 
conditions show a very strong, near surface stratification occurring within 
the water column, ideal for internal tide generation and evolution of 

solitons.16 

Figure 10. CTD data taken from NAVO sponsored Calcofi cruises. Location of cast: 

36.74°N, 122.027°W. (a) February 2001. Weak, deep stratification is 
clearly observed in the Temperature and Density Anomaly vertical 
profiles. The mixed layer can be observed extending down to 50m depth. 

N vs. depth shows Nmax occurring much at 50m with a magnitude of 0.01 
s' 1 , (b) January 2002. The Temperature and Density Anomaly vertical 
profiles identify an almost non-existent mixed layer at approximately 60m. 
Again, Nmax defines the mixed layer at 60m with a magnitude of 0.02 s' 1 . 
Winter conditions show a very week, deep stratification in the water 
column, not ideal for the generation of an internal tide and solitons.17 
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(a) Ebb flow across a shelf break causes a depression in the pycnocline. 

(b) Formation of a steep edge shoreward bore during slack tide, (c) 
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edge continues to steepen due to dispersive effects, (d) The steep leading 
edge of the bore degenerates into solitons through the dispersive properties 

of internal waves.19 
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Location of Monterey Bay and the MISO platform off Monterey, CA. The 
star represents the location of the MISO array. The axis indicates the CS 
normal coordinate system used in the MISO bottom boundary layer study, 

with +x offshore, -x onshore and +y south AS and -y north AS.28 

Yearday 228-228.15. (a) Temperature profile timeseries. (b) ADCP CS 
velocity profile timeseries. (c) Current vectors. Vertical alignment of 
isotherm and current displacement are indication of a soliton and internal 

tidal bore.34 

YD 229-229.1 (a) Temperature profile timeseries. (b) CS ADCP current 
velocity profile timeseries. (c) Current vector timeseries. CS ADCP 

current velocities indicate solitons as part of a wave train.36 

YD 228-229. (a) Temperature profile timeseries. (b) CS ADCP velocity 
profile timeseries (c) Low pass filtered current vectors. Although the 
temperature timeseries has missing data, the semidiurnal internal tidal bore 

is evident in the CS ADCP velocity and current vector timeseries.39 

YD 229-230(a) Temperature profile timeseries. (b) CS ADCP velocity 
profile timeseries (c) Low pass filtered current vectors. Although the 
temperature timeseries has missing data, the semidiurnal internal tidal bore 
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known as spring tides and period of minimum amplitude are neap tides. 
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I. INTRODUCTION 


Surface disturbances caused by strong, nonlinear, near-surface internal waves are 
commonly seen over the continental shelf as bands of slick and rough water propagating 
shoreward (Kropfli et al 1998). These internal waves, or solitons, evolve on the leading 
edge of an internal tidal bore that has steepened due to dispersion effects. As these 
solitons propagate shoreward and begin to shoal, their form changes, and they start to 
affect the bottom boundary layer above the ocean bed. From their first observation in 
1834 by John Scott Russel as a surface disturbance propagating down a canal, few studies 
have been conducted on their effect on the bottom boundary layer. Past studies have 
primarily been focused offshore and on their effects on the current and temperature fields 
as well as their dynamics. Studies conducted by Baines (1972) in the Mid-Atlantic 
Ridge, Apel (1985) in the Sulu Sea, and Holloway (1987) off northwest Australia were 
mainly focused on internal tide generation over a shelf break and on characterizing the 
resulting solitons. More recent studies conducted by Liu (1988) and Chapman et al 
(1997) in the New York Bight and Stanton and Ostrovsky (1998) off northern Oregon 
began to focus on solitons as they propagated away from their generation site. Still, little 
has been discussed on soliton effects further inshore on the shelf and their effects on the 
bottom boundary layer due to shoaling. Liu (1988) briefly touched on the shoaling 
parameter as it relates to soliton growth and decay. The purpose of the inner shelf soliton 
experiment conducted in Monterey Bay during the summer of 2002 was to closely 
examine soliton and internal tidal bore effects on the bottom boundary layer due to 
shoaling at an inner shelf site only 600m from the surf zone. 

Certain environmental conditions are necessary to allow an internal tide and 
solitons to develop in coastal regions. Strong, shallow stratification in the water column, 
which frequently approximates a two-layer system, provides the background for these 
waves to form (Stanton and Ostrovsky 1998). A significant cross-shelf barotropic tidal 
current provides the necessary forcing in order to generate an internal tide. This internal 
tide propagates shoreward as an internal tidal bore, and the baroclinic structure of the 
bore causes water near the surface to be pushed onshore and water nearer the bed to be 
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pulled offshore. As the leading edge of the tidal bore steepens it begins to degenerate 
into solitons due to the dispersive properties of internal waves. These solitons have the 
same baroclinic structure as the tidal bore, but with much shorter wavelengths, and cause 
currents to flow in the same pattern. The main objective in the bottom boundary analysis 
is to prove that the strong push and pull of water associated with tidal bores and solitons 
are felt near the surface and all the way down the water column to the bed. From this it 
can be shown if stresses generated by offshore movement of water near the bed act 
together with near bed forcing from surface gravity waves to disrupt the bed and cause 
sediment to be suspended and thus, produce a significant net transport. 

The location and complexity of the shelf break plays a significant role in the 
propagation of internal bores. As in previous studies conducted, the location of the shelf 
break was approximately parallel to the coast causing Stanton (1998), Liu (1988), 
Chapman (1997) and others to conclude that tidal bore and soliton propagation are 
largely cross-shelf, with currents moving on and offshore. However, due to the 
complexity of the shelf break and canyon within Monterey Bay, it is likely that internal 
tidal bores might have significant alongshore (AS) components. It is hypothesized that 
the morphology of the Monterey Canyon to the north of our observation site has caused 
this change in direction and has funneled the internal tide up through the canyon axis and 
refracted it along its steep walls. 

The Monterey Inner Shelf Observatory (MISO) was designed and implemented 
by Stanton (1999) to provide a long-term observation capability to study the evolution of 
a sandy bed at an inner shelf site. The objective of this paper is to analyze the effects of 
tidal bores and solitons on the bottom boundary layer by measuring stress components, 
sediment suspension profiles, net transport, and to characterize the solitary internal wave 
and internal tidal bore summer climatology at the site. 

A. BACKGROUND 

A regional study on the offshore and onshore oceanographic and climatological 
conditions of Monterey Bay was conducted for the year prior to bottom boundary layer 
observations. The main objective of this regional study was to look at seasonal scale 
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changes in the oceanographic structure that determine the presence or absence of internal 
tidal bores within the bay. 

1. Monterey Bay 

Monterey Bay is located along the central California coast between 121°W- 
123°W longitude and 35°N-37°N latitude (Figure 1). The bathymetry offshore of central 
California is complex, with canyons and ridges dominating the underwater topography as 
seen in Figure 1. MISO is located in the southern end of Monterey Bay on the inner 
shelf. North of the MISO site and at the center of the bay is the Monterey Canyon, which 
reaches depths of up to 4000m. The Monterey Canyon is a classic example of a 
submarine canyon. Several studies have been conducted to determine the effects on the 
internal tide within the canyon (Petruncio et al 1997). As the canyon winds its way into 
Monterey Bay, it cuts into the shelf break with its steep walls and complex bathymetry. 
The southern edge of the Monterey Canyon wall provides a steep shelf break area north 
of the MISO site, with depths jumping from 1000m to 100m over a short 



Figure 1. Location of Monterey Inner Shelf Observatory (MISO) site. The Monterey 
Canyon can be seen bisecting the bay with its steep walls and complex bathymetry. The 
‘x’ denotes the location of MISO and the 3 bars indicate possible generation sites for 
internal bores. The 2 stars indicate buoys Ml (inshore) and M2 (offshore). Courtesy of 
MBARI. 
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distance of 5000m. Roughly linear segments of rapid off-shelf depth changes are 
indicated in Figure 1 with the shelf break directly offshore approximately 10km 
northwest of MISO. These three shelf break areas can be considered as potential 
generation sites for internal tidal bores within the bay. 

2. Offshore Ocean Observations 

Data collected by the Monterey Bay Aquarium Research Institute (MBARI) was 
used to characterize the offshore stratification of the bay for the year (2001-2002) 
previous to the July and August 2002 bottom boundary layer observation made at the 
MISO site. Temperature, salinity, and pressure from in-situ sensors separated by 20m on 
the Ml (inshore 36.75°N, 122.03°W) and M2 (offshore 36.70°N, 122.39°W) moorings 
(Figure 1) were analyzed and used to calculate potential density (c 0 ). The data were 
collected at 20 sec intervals and averaged to 2-day samples in this analysis. From this, 
the potential density was then used to calculate N, the Brunt Vaisala or buoyancy 
frequency, 

/v=[-Ap (s' 1 ) (1) 

P Clz 

d p /dz was calculated using the method of finite differencing. MBARI salinity, 
temperature, and pressure data were collected every 20m down to depths of 200m (Ml) 
and 300m (M2) at 20sec intervals. These data were linearly interpolated to depths of 
10 m and used to estimate d p /dz using finite differencing. 

dp = P;-P ; -1 (2) 

dz Z z -Z z _ 1 

where subscript z denotes the depth level and subscript zl denotes the level above z. 
This method, however, has certain limitations. Unfortunately, the 20m sensor spacing of 
Ml and M2 under samples sharp density gradients with scales « 20m. Therefore, this 
estimate of N can be regarded as a lower bound. More significant is that the mixed layer 
can nearly reach the surface, so the 20m depth spacing would not be able to resolve this 
shallow stratification. Due to this limited vertical sampling, the N calculated here is a 
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smooth, coarse estimate of the actual stratification; however it does give a long-term 
timeseries of coarse upper ocean structure. 

Together, o 9 and N are used to evaluate the seasonal stratification cycle of the 
bay and to determine if conditions can support the generation of an internal tide and 
solitary internal waves. 

3. Winds 

Climatological wind data measured at buoys Ml and M2 show hat winds are 
predominately from the northwest at 5-10 ms" 1 during most part of the year in Monterey 
Bay. There are brief interruptions in the daily wind pattern during the winter months due 
to winter storms moving through the area (Figure 2). The summer period is strongly 
dominated by persistent northwest winds. Daily wind patterns drive the surface 
circulation of the bay and have a large effect on surface layer mixing and upwelling 
sources. Near surface stratification can rapidly change due to a reversal of the winds. 
When winds are from the north, they transport with them cold, fresh water, in contrast to 
when they are from the south, when they transport warm, more saline water up towards 
the bay. Depending on the wind patterns, and direction of the surface currents, the 
presence of shallow stratification can quickly change. So consequently, winds may 
impact the formation and propagation of an internal tide and solitons. 

Wind forcing during the spring and summer months also plays a significant part 
in coastal upwelling and the advection of water into and out of the bay (Rosenfeld el al. 
1993). This too can have a strong impact on the stratification of the bay. Collins et al. 
(2002) noted that surface waters in Monterey Bay were coolest during the spring/summer 
upwelling period. He also noted coldest temperatures (<10°) at 80 dbar during periods of 
intense upwelling. In contrast, the winter months typically see a decrease in winds from 
the northwest. Along the coast, upwelling begins to cease and lie Davidson current 
begins to surface (Pennington and Chavez 1999). Winter storms move through the 


5 



1 o 


Ml WtMtJ 



y«»rt)»y 200* 


Ml buoy Nov2001 -Oct20D2 



spring-► IM^- 

vvrrt pf-» summer—» 


M2 buoy Ap*20D1 •May2Q02 



«prmg ► faff—»■ spring ► 


summer^ 


winter ► 




Figure 2. Wind speed and direction from MBARI Ml and M2 buoys. Data is averaged 
every 7 days over the period of 1 year. Scalar averaging makes winter storms much more 
apparent as seen by the arrows. During the winter period, these winter storms disrupt the 
normal wind patterns and help to lower the mixed layer in the water column. They help 
to mix the upper part of the water column and decrease the insulation. These factors 
combine to deepen the stratification during the winter. 
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region, bringing with them intense winds that result in a deepening of the mixed layer. 
This turbulence decreases the insulating effect in the upper part of the water column. 
Cessation of upwelling and the start of winter storms combine to decrease the near 
surface stratification in the bay during the winter months. Collins (2002) identified 
surface waters to be warmest in the fall and on into the winter. He observed waters below 
80 dbar to be warmer during the winter with some temperatures exceeding 14°C. The 
stratification maximum is typically found deeper in the water column and relatively 
weaker during the winter period. 

4. Circulation 

Wind driven currents produce a strong surface circulation within and just offshore 
of Monterey Bay. Since winds are predominantly from the northwest in the summer, the 
surface current flows from north to south as it passes the entrance to Monterey Bay. 
During upwelling conditions, the source of this water is located just north of the bay 
centered at Pt. Ano Nuevo (Rosenfeld el al. 1993). This cold, salty water typically 
dominates the water mass source in the bay and helps to sustain the strong near surface 
stratification that is present in the summer. Strub et al. and Chelton el al. (1987, 1988) 
have identified a seasonal cycle associated with currents in Monterey Bay. They 
observed a southward flow during the upwelling season (spring, summer) and a 
northward flow the rest of the year. During the spring/summer upwelling period the 
water is cooler at 80 dbar due to a southerly flow in contrast to the winter when water is 
warmer at the same depth due to a northerly flow. 

Within the bay, as winds carry surface water down from the north, surface 
circulation is observed to move in a counterclockwise direction, advecting the cold, salty 
water into the bay. Just outside of the bay, the mean surface circulation is in a clockwise 
direction. 

5. Stratification 

Stratification during the summer is stronger than in the winter within Monterey 
Bay as seen daily in sampled temperature and density profile timeseries (Figures 3 and 
4). MBARI data show strong, near surface stratification present offshore at Ml and M2 
moorings during the summer months. In the summer there are strong near-surface 
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re 3. (a) Temperature timeseries of Ml mooring in Monterey Bay. 13Nov01-170ct02. 

(b) Salinity timeseries of Ml mooring in Monterey Bay. (c) Sigma-theta (Oe) timeseries 
of Ml mooring in Monterey Bay. Strong, shallow stratification can be seen during the 
summer months. The winter months show a weaker and deeper stratification. Conditions 
during the summer months are more ideal for the generation of an internal tide and 
solitions. 
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:e 4. (a) Temperature timeseries of M2 mooring in Monterey Bay. 12Apr01-08May02. 

(b) Salinity timeseries of M2 mooring in Monterey Bay. (c) Sigma-theta (Oe) timeseries 
of M2 mooring in Monterey Bay. Strong, shallow stratification can be seen during the 
summer months. The winter months show a weaker and deeper stratification. Conditions 
during the summer months are more ideal for the generation of an internal tide and 
solitions. 








thermal gradients and the stratification nearly reaches the surface (Figure 5). At M2, 
strong near surface stratification is observed as well (Figure 6). In each case, thermal 
gradients extend above 40m depth. Distinct undulations in the warm upper layer can be 
observed, particularly in the M2 data. This may be associated with changes in forcing 
during the spring and summer in which the water column restratifies every 2-3 weeks due 
to a reversal of the winds which intermpts the upwelling conditions. High temperatures 
and low salinity observed at the surface both contribute to a strong near surface 
stratification. Profile timeseries of salinity and potential density for Ml and M2 are 
similar in structure to the thermal structures in Figures 5 and 6. Figures 5 and 6 clearly 
show the strong density gradients that are present near the surface during the summer. 
The summer months distinctly show the presence of a strongly stratified water column 
near the surface, which approximates a two-layer system. However, due to the 20m 
sensor spacing on the Ml and M2 moorings, the density structure is unable to be resolved 
to finer scales. This limits our estimate of the Viasalla frequency to a lower bound, and 
prevents a definition of mixed layer depth. 

During the winter months, however, near-surface stratification disappears almost 
completely (Figures 7 and 8). At both the Ml and M2 moorings, stratification is weaker 
and located much deeper in the water column. The deeper mixed layer during the winter 
is due to an increase in surface stress from winter storms and the cessation of upwelling. 
Despite having low average winds, strong, intermittent winter storms help to deepen the 
mixed layer with sudden increases in colder winds. Because Figure 2a is a vector 
average over 1 year, winter storms would have a weak signature, as they are episodic, 
multi-directional events. However, when scalar wind speeds are averaged, the winter 
storms become much more apparent (Figure 2). The winter months also show the 
absence of a near-surface salinity stratification, indicating that the mixed layer is deeper 
in the water column as seen in Figures 7 and 8. 

Profile timeseries of N (Brunt-Visailla frequency) vs. depth help to determine the 
strength of the water column stratification. Strong stratification near the surface, seen in 
the summer average (Figures 5 and 6), produced a lower bound N max value of 0.014 s' 1 at 
a depth of 50-60m. Weaker stratification, which is seen in the winter, yields values of 
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Figure 5. Summer conditions at Ml mooring in Monterey Bay for yearday 172 to 266. (a) 
Temperature profile indicates a strong and shallow stratification in the water column, (b) 
Salinity profile shows a strong and shallow stratification as well, (c) Sigma-theta (Oe) 
profile agrees well with the temperature and salinity profiles showing a strong and 
shallow stratification, (d) N (buoyancy frequency) calculated as a summer mean using 
20m Oe differences gives indication that shallow stratification is occurring in the water 
column. The strong, shallow stratification in the summer is ideal for the generation of an 
internal tide and solitons. 
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Figure 6. Summer conditions at M2 mooring in Monterey Bay for yearday 172 to 266. (a) 
Temperature profile indicates a shallow and strong stratification in the water column, (b) 
Salinity profile shows a strong and shallow stratification as well, (c) Sigma-theta (Oe) 
profile agrees well with the temperature and salinity profiles showing a stronger and 
shallower stratification, (d) N (buoyancy frequency) calculated as a summer mean using 
20m Oe differences and plotted vs. depth confirms that the water column has a very 
strong and shallow stratification. The strong, shallow stratification in the summer is ideal 
for the generation of an internal tide and solitons. 
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MBARI Ml winter 2001-2002 




Figure 7. Winter conditions at Ml mooring in Monterey Bay for 2001-2002 yearday 356 to 
80. (a) Temperature profile indicates a very deep and weak stratification in the water 
column, (b) Salinity profile shows a deeper stratification as well, (c) Sigma-theta (Oe) 
profile agrees well with the temperature and salinity profiles showing a deep and weak 
stratification, (d) N (buoyancy frequency) calculated as a winter mean using 20m Oe 
differences plotted vs. depth clearly indicates that stratification is weak and occurs very 
deep within the water column. The weak, deep stratification in the winter is not ideal for 
the generation of an internal tide and solitons. 
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Figure 8. Winter conditions at M2 mooring in Monterey Bay for 2001-2002 yearday 356 to 
80. (a) Temperature profile indicates a very weak and deep stratification in the water 
column, (b) Salinity profile shows deeper stratification as well (c) Sigma-theta (Oe) 
profile agrees with the temperature and salinity profiles indicating that the water column 
has a weak and deep stratification, (c) N (buoyancy frequency) calculated as a winter 
mean using 20m Ge differences plotted vs. depth clearly indicates that stratification is 
present, but very deep in the column. The weak, deep stratification in the winter is not 
ideal for the generation of an internal tide and solitons. 
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Nmax on the order of 0.008 s’ 1 , significantly less than those of summer months. These 
values of N max are generally found deeper in the water column at a depth of 80m (Figures 
7 and 8). 

Together these timeseries give a clear indication that the water column during the 
summer has a strong, near surface stratification. This makes conditions during the 
summer more ideal for the generation of an internal tide and solitons. Because the 
MBARI data have limited resolution of the water column due to its 20m sensor spacing, 
CTD data were taken from NAVO sponsored Calcofi cruises within Monterey Bay in 
order to validate and support our conclusions on the stratification. The 1000m CTD casts 
from July 2000 and August 2001 show a finer resolution of the strong near surface 
stratification. Conversely casts from February 2001 and January 2002 near the mooring 
sites show deeper, weaker stratificaiton. These CTD casts provide a 2m resolution of the 
temperature, salinity, and density structures. The CTD casts were made in the vicinity of 
the Ml mooring at 36.73°-36.74°N and 122.028°-122.026°W. Figures 9 and 10 show the 
CTD casts in which vertical profiles of temperature, salinity, density anomaly, and N 
were computed. Figure 9 depicts the summer conditions, showing a mixed layer at l-2m 
depth for July 2000 and at 10m depth for August 2001, whereas Figure 10 depicts the 
winter conditions with a mixed layer at 50m for February 2001 and January 2002. When 
compared with N, we can see better resolved values of the stratification, with N max for the 
summer of 0.03 s’ 1 near the surface at l-2m depth. This value is considerably larger than 
Nmax for the winter which had a value of 0.015 s’ 1 at a depth of 60m. 

6. Tides 

A barotropic semidiurnal tide is present in Monterey Bay, resulting in two highs 
and two lows that occur daily. The tides are linked to astronomical forcing dominated by 
the position of the sun and moon relative to the earth. During the summer, the new and 
full moons generally fall around the middle and end of each month, with each month 
changing slightly due to the lunar cycle and sidereal day versus the calendar day. 
respectively. It is at these times when the earth, sun, and moon are aligned that the tides 
are very large and are known as spring tides. When the moon is in first or last quarter 
(beginning and third week of each month), the tides are generally weaker and are known 

as neap tides. These tides occur when the earth, sun, and moon are in quadrature. Spring 
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Figure 9. CTD data taken from NAVO sponsored Calcofi cruises. Location of cast: 
36.74°N, 122.027°W. (a) July 2000. Strong, shallow stratification is clearly observed in 
the Temperature and Density Anomaly vertical profiles. The mixed layer can be 
observed at l-2m depth. N vs. depth shows Nmax occurring at the surface (l-2m) with a 
magnitude of 0.03 s' 1 , (b) August 2001. During late summer, stratification starts to move 
deeper into the water column. It can be found at 10m depth. Nmax is also seen to 
decrease to 0.015 s 1 and move down as well (10m). Summer conditions show a very 
strong, near surface stratification occurring within the water column, ideal for internal 
tide generation and evolution of solitons. 
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Figure 10. CTD data taken from NAVO sponsored Calcofi cruises. Location of cast: 
36.74°N, 122.027°W. (a) February 2001. Weak, deep stratification is clearly observed in 
the Temperature and Density Anomaly vertical profiles. The mixed layer can be 
observed extending down to 50m depth. N vs. depth shows Nmax occurring much at 
50m with a magnitude of 0.01 s' 1 , (b) January 2002. The Temperature and Density 
Anomaly vertical profiles identify an almost non-existent mixed layer at approximately 
60m. Again, Nmax defines the mixed layer at 60m with a magnitude of 0.02 s' 1 . Winter 
conditions show a very week, deep stratification in the water column, not ideal for the 
generation of an internal tide and solitons. 
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tides have offshore tidal surface currents up to 20 cms" 1 , while neap tides have tidal 
surface currents up to 10 cms" 1 , as observed by Petruncio (1993) and Paduan et al. (1995) 
using CODAR. They observed tidal surface currents to be larger at the head of the 
canyon (15 cms" 1 ) and to be aligned with local bathymetry at the southern end of the bay. 
The K1 and 01, which are the lunisolar diurnal constituents that expresses the moon’s 
declination, have periods of approximately 24 hours (NOAA). The M2 and S2, which are 
the principal lunar and solar semidiurnal constituents, have periods of approximately 12 
hours. Together, each of these tidal constituents contributes to over 80% of the total 
amplitude and tidal forcing in Monterey Bay. 

A barotropic semidiurnal tide dominates the bay and is characterized as a surface 
wave with particle motion following the high and low tides. Theses tides have extremely 
long wavelengths on the order of thousands of kilometers. The astronomical forcing of 
M2 was observed to be the strongest, contributing the most to the barotropic tidal activity 
within the bay (Petruncio et al. 1997). The second strongest constituent observed was the 
S2, which was 26% of M2. 

7. Internal Tide Generation Mechanisms 

With a strong, shallow stratification present, a sharp shelf break 10km offshore, 
and a strong barotropic cross-shelf current, conditions in Monterey Bay are ready for the 
generation of an internal tide. The internal tide is generated in the following manner, 
following Chapman et al (1997): 

As the ebb tide begins to flow out of the bay and makes contact with the shelf 
break, it causes a depression in the isopycnals (Figure 11a). This isopycnal depression, 
or thickening of the upper layer of the stratified column, is the result of a standing lee 
wave formed by offshore flow over the hydraulic jump at the shelf break. When the flow 
ceases during slack tide, the pycnocline depression is able to propagate internal tidal 
bores both onshore and offshore (Figure lib). During the flood tide, the steep edge 
onshore bore propagates shoreward (Figure lie) (Holloway 1987). The leading edge of 
this internal bore will begin to steepen due to dispersion effects as the bore continues to 
propagate shoreward. Nonlinearities will begin to arise as the leading edge continues to 
become steeper. These nonlinearities are what give rise to solitons (Figure lid) (Liu 
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Figure 11. (a) Ebb flow across a shelf break causes a depression in the pycnocline. (b) 

Formation of a steep edge shoreward bore during slack tide, (c) Internal bore propagates 
shoreward assisted by the flood tide. Leading edge continues to steepen due to dispersive 
effects, (d) The steep leading edge of the bore degenerates into solitons through the 
dispersive properties of internal waves. 
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1988). As the internal tidal bore continues propagating shoreward across the shelf, the 
solitons begin to grow in size and energy. 

The generation of an internal tide is closely related to the slope of the shelf break 
and tidal forcing of the and S 2 frequencies (Baines 1982). The rate of conversion of 
surface tidal energy into internal tidal energy is most notable over continental slopes and 
increases rapidly with the topographic height (Baines 1972). Baines (1985) modeled 
internal tide generation in laboratory experiments as a function of the amplitude of the 

u , 

barotropic tidal forcing given by the Froude Number, F = —^—, where u sb is the 

c 

XV 

maximum barotropic tidal velocity at the shelf break and c M ,is the long wave speed of the 
lowest internal mode. He found a strong correlation between the slope of the shelf and 
the Froude number. For a two-layer stratification, large Froude numbers (>1) were 
observed indicating that linear wave theory was not adequate to predict internal tide 
generation over the hydraulic jump that occurs seaward of the shelf break during the ebb 
phase of the barotropic tide. However, in his laboratory experiments, he was able to 
conclude that for high Froude numbers and steep slope geometry, a hydraulic jump will 
form seaward of the shelf break during the ebb phase of the tidal cycle. For an 
approximate 2 -layer stratification with one layer thicker than the other, an internal bore 
will form. Values of F>\ are commonly seen in the open ocean over continental shelves 
and help to explain the nonlinear character of internal waves of tidal origin (Baines 
1985). 

Holloway (1999) also came to the conclusion that internal tide generation is 
strongly dependent on the relationship between the slope of the bathymetry and the slope 
of the internal wave characteristics. He defined the slope characteristics (s) as the path 

. v • v v + Ico 2 -/ 2 , . , 

along which the wave energy propagates as s = — ,—-- 7 , where (0 is the wave 

-VAV- o)" 

frequency, / is the Coriolis parameter, and N is buoyancy frequency. He further defined 

oc oc 

critical values for bathymetry slope (a ): critical when —= 1 , subcritical when —< 1 , 

s s 
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and supercritical when —>1. He observed that little generation occurs for weak 

s 

subcritical slope, while resonance near the sea bed may occur for critical slopes. 

The continental shelf break 10km directly offshore of the MISO site, as well as 
the southern wall of the Monterey Canyon to the north, are shelf breaks for the generation 
of an internal tide within Monterey Bay. Because the continental shelf break is directly 
offshore and parallel to the coast, it will produce internal bores that propagate in a cross¬ 
shore (CS) direction, whereas internal tidal bores generated along the southern wall of the 
canyon would likely produce more alongshore (AS) directed bores that could be refracted 
by the shallow, inner shelf topography. 

B. SOLITONS 

Solitons move in packets frequently bound to the steep front of an internal tidal 
bore. They are widely observed in coastal regions when the water column is stratified and 
generally propagate CS and can be found to be moving on and offshore (Apel 1983). 
They are strongly nonlinear internal waves that are coupled with the internal tide and 
generate strong current pulses in the water column. 

Solitons can be characterized by their large amplitudes, vertical displacements of 
the isotherms, and strong currents associated with them. They are commonly observed as 
part of a wave train on the leading edge of the steepening internal bore and are readily 
identified by the vertical alignment of current and isotherm displacements (Liu 1988). In 
a very strongly stratified coastal region, Stanton and Ostrovsky (1998) observed vertical 
displacement amplitudes of a soliton to be 2-4 times greater than the pycnocline depth 
with lateral currents up to 1ms" 1 off the coast of Northern Oregon. Within these wave 
packets, there are rapid, unidirectional pulses in the direction of propagation of the 
solitons. In a coastal region with much weaker stratification, Chapman et al (1997) 
observed typical speeds of l-2cms _1 , with some peaks as great as 5cms _1 . They also noted 
large temperature variations of up to 6-7°C within the New York Bight. In both these 
regions a distinct shelf break was located offshore of the solitary internal wave 
observation site, similar to the topographic structure in the southern part of Monterey 
Bay. 
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When the leading edge of the internal bore steepens due to dispersion effects, the 
amplitude of the solitons increases as well. As the amplitude increases, the width of the 
soliton decreases (Stanton and Ostrovsky 1998). Increasing order of nonlinear functions 
are then needed to describe their structure. Stanton and Ostrovsky (1998) defined the 
width of a soliton as the time interval between displacements at half the maximum 
amplitude. Solitons off northern Oregon were observed to have durations of 180s. For 
the 25 solitons observed, Stanton and Ostrovsky (1998) noted that as the amplitudes 
varied from 25m to 5m, the durations (and here widths) remain generally unchanged at 
180s. 

Soliton displacement can occur near the surface or near the bottom of the water 
column depending on the background stratification, which is itself modulated by the 
internal tide on the inner shelf. As a near surface soliton propagates shoreward and 
begins to shoal, its form changes and may affect the bottom boundary. Liu (1988) notes 
that as the wave train propagates inshore from the shelf break at a depth of 90m to the 
nearshore observation site at a depth of 30m, the linear wave speed is reduced to one-half 
its initial value, dispersive effects decrease as the wave is compressed, nonlinear effects 
decrease, and shoaling increases. Liu (1988) also observed that a critical depth existed at 
45m depth in the NY Bight at which point the wave packet propagates onto the shelf and 
disappears. From this, he concluded that varying depth on the evolution of solitons fell 
into two categories: 1) The linear wave speed is modified along the path of wave 
propagation that leads to spreading out or compression of the waves that disintegrate into 
a dispersive wave train at the critical depth. 2) The slope of the shelf produced a shoaling 
term proportional to the wave amplitude. This shoaling was observed to increase the 
soliton’s amplitude and increase the wave energy. However, Liu (1988) also noted that 
turbulent processes begin to erode the peaks of the well-developed soliton, decreasing 
their amplitude, and increasing their width. The radial spreading and dissipation effects 
eventually suppressed the shoaling effects and it was found that the evolution of solitons 
was a balance between nonlinear and dispersive effects. The eddy viscosity present 
eroded the sharp peak of the lead soliton causing it to decay, while the shoaling effect 
balanced the dissipation effect for the rest of the solitons in the wave packet allowing 
them to grow. 
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D. FIELD EXPERIMENT 


A. PHYSICAL SETUP 

The field data used in this study are collected at the Monterey Inner Shelf 
Observatory (MISO) site, which is part of the Rapid Environment Assessment Laboratory 
(REAL), located off of Del Monte Beach and the Naval Postgraduate School (NPS) in the 
southern end of Monterey Bay. MISO was designed and implemented by Tim Stanton 
and his research group at NPS. It is sponsored by the Oceanography Department at the 
Naval Postgraduate School and by the Office of Naval Research as part of the Coastal 
Processes and Physical Oceanography programs and the Shoaling Waves Experiment. 
The facility was deployed in January 2000 and is ongoing. The inner shelf soliton 
experiment was a two-month study that took place during the 2002 summer, with rapid 
CTD profiles supplementing the detailed current measurements for 2 weeks within the 
period. 

The MISO instrument frame is located 600m offshore from the beach in 12m of 
water and lm above the bed (Figure 1). Long timeseries of high resolution 
measurements made on the MISO frame are possible because of the high bandwidth data 
communications and power links between MISO and the laboratory onshore. These 
instruments are being used to study the interaction of winds, waves, and sediment bed on 
the inner continental shelf, just offshore from the surf zone. A detailed diagram of the 
frame and instruments is shown in Figure 12. 

B. INSTRUMENTS 

1. Paroscientific Pressure Transducer 

A high precision digital pressure sensor (Paro-scientific) measures the height of 
water above the frame by sensing the hydrostatic pressure of the water column and 
propagating waves to produce a record of tides and surface gravity waves. It is located 
approximately 1.2m above the bed. In this study, 5 minute averages of pressure are used 
to filter out surface gravity waves. The pressure sensing element in this instrument is an 
oscillating quartz crystal that provides very high long term stability and resolution. 
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Figure 12. MISO instrument frame as viewed from offshore looking onshore. The 
instruments include (from left to right): camera component of the Structured Light and 
Camera (SLAC), Bistatic Coherent Doppler Velocity and Sediment Profiler (BCDVSP), 
structured light component of the SLAC, Scanning Acoustic Altimeter (SAA), Paro- 
scientific pressure sensor, and Braodband Acoustic Doppler Current Profiler (BADCP). 


2. ADCP 

An upward looking WorkHorse Acoustic Doppler Current Profiler (ADCP) 
manufactured by RD Instruments measures the water current every 0.5m over 27 equally 
spaced bins from the frame to the near ocean surface. The ADCP starts each ping cycle 
by transmitting a broadband acoustic pulse into the water column along each of the four 
narrow acoustic beams that are tilted 20° off the vertical. Each transducer detects 
Doppler shifts of scatterers in the water column moving with a component of velocity 
aligned with the beam direction. These Doppler shifts are estimated every 0.5m in the 
MISO setup by range-gating the signals. The resulting set of four radial Doppler shifts at 
each range bin are then resolved into three components (u, v, w) of velocity assuming 
there is no variation in velocities across the diverging beams. This limits the accuracy for 
wind wave-induced velocities at the most distant bins near the surface; however, in this 
study lower frequency motions are being considered. Acoustic reflection caused at the 
surface frequently results in contamination of the upper 1-2 bins. The ADCP is also 


24 










equipped with a heading and tilt sensor allowing the instrument coordinates to be 
transformed into earth-leveled. North referenced coordinates. 

3. BCDV 

A high resolution Bistatic Coherent Doppler Velocity and Sediment Profiler 
(BCDVSP) developed at NPS by the ocean turbulence research group measures the three 
component velocity vector and sediment concentration every 0.01m over 100 equally 
spaced bins from a range of 0.60m above the bed (Stanton 1999). The BCDV transmits a 
1.25MHz acoustic signal, and uses coherent sampling and a bistatic geometry to receive 
backscattered energy. The BCDV backscatter energy is proportional to the sediment 
concentration, over a range of from 0.01 to lOOg/L. These small- scale measurements are 
used to define the mean and turbulent structure of the water column above the sandy bed, 
which results from currents interacting with the bed. 

4. CTD 

A Direct Reading Micro CTD developed by Falmouth Scientific Inc. was used to 
collect temperature, salinity, and pressure data over 127 bins equally spaced every 0.1m. 
Conductivity was measured through a FSI inductive cell with an accuracy of +/- 0.002 
S/m and temperature was measured with a platinum resistance thermometer with an 
accuracy of +/- 0.002°C. Pressure was measured through a micro-machined silicon 
pressure sensor with an accuracy of +/- 0.02%. An automated, computer controlled 
winch developed at NPS repeatedly profiled the CTD from 2m HAB to the surface every 
2 minutes. The CTD provided temperature and salinity profile timeseries to measure the 
evolution of the stratification of the water column. Poor temporal response of the 
conductivity cell severely limited the vertical resolution of the conductivity, and hence, 
salinity profiles during these observations. Consequently, isotherm displacements are 
used instead of isopycnals in the following analysis. 

5. Meteorology 

Wind speed and direction from the Marine Operations Lab MET instrument tower 
were used as a measure of the local winds. The 8m tower was located approximately 
700m directly shoreward of the MISO array. 
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HI. DATA PROCESSING 


A. COORDINATE TRANSFORMATION 

In order to correctly analyze the field data collected by these instruments, the 
coordinate system of the ADCP and BCDV were resolved to a cross-shore (CS) normal 
coordinate system with +x offshore and +y along-coast toward Monterey (see map. 
Figure 13). The ADCP is oriented such that it is in magnetic coordinates where +/-u is 
east/west and +/-v is north/south. The ADCP heading (ADCPhdg 0 ) measured by the 
internal compass was 309°-312° magnetic. Shore normal to the beach was estimated to 
be 320° magnetic. To provide a coordinate system with +x offshore, -x onshore, and +y 
alongshore (AS) towards Monterey, and -y alongshore (AS) towards Marina, the ADCP 
heading must be rotated an additional -90° (Figure 1) 


ADCP rotation 0 = ADCP h dg° + (320° - ADCP hdg 0 ) - 90° (3) 


The BCDV is oriented in a frame coordinate system, 218° from magnetic North, 
with a heading of 308°magnetic. A 12° rotation had to be applied to resolve it to a shore 
parallel coordinate system. To further resolve it to a CS normal coordinate system, a -90° 
was applied to the data. 


BCDV rotation 0 = 12° -90° 


( 4 ) 


To analyze the effects of wind forcing, wind data had to be resolved into the CS 
normal system as well. Wind from the onshore MET tower was measured in a north- 
south coordinate system in degrees True. Because the wind was measured in a true 
coordinate system a 15° correction had to be made for the magnetic variation. Wind 
speed and direction were used to calculate the u and v components of the wind. These 
components were then transformed into the CS normal coordinate 
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Figure 13. Location of Monterey Bay and the MISO platform off Monterey, CA. The star 
represents the location of the MISO array. The axis indicates the CS normal coordinate 
system used in the MISO bottom boundary layer study, with +x offshore, -x onshore and 
+y south AS and -y north AS. 
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system by applying a heading rotation to the u and v components to find the CS (windcs) 
and long-shore (windLs) components. 


wind cs = windspeed (—) * cos[ winddirect ion - 15" ( rad) - 180 * ] (ms -1 ) 


pi 


180 


(5) 


wind^ = windspeed (—) * sin[ winddirection — 15 0 ( rad) - 180 * I ( ms ') 


( 6 ) 


wind rotation =360 o +(320°-360 o >90° 


(7) 


B. HEIGHT ABOVE BED REFERENCE 
1. CTD 

Once the coordinate system was resolved, a fixed vertical depth reference had to 
be established prior to analyzing any of the field data collected by the CTD and velocity 
profiling instruments. As the focus of this study is the bottom boundary layer, height 
above the bed (HAB) reference in meters was selected as vertical coordinate. The 
pressure-referenced CTD profiles were converted to this vertical reference by using the 
MISO pressure sensor. The height in meters above the bed of the paroscientific was a 
known value of 1 . 2 m and the pressure of the paroscientific (p paro ) was taken to be the 
pressure measured by the instrument (p pa mtatoi) minus the barometric pressure (pbarometric) 

Pparo ~ P parototal Pbarometric (dbar) ( 8 ) 


Using the hydrostatic equation and converting the 5 minute pressure averages from dbar 
to Pa, the height above the paroscientific (z pa m) was then calculated. 
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(9) 


kg ni 

P = p(-y)g(— )z(m) (dbar) 

m s 


= P„ 


10 4 . . 

*- (m) 

P£ 


( 10 ) 


The distance below the surface (zq) to the CTD was then calculated in the same manner. 


7 = 


_ WF CTD 


Pba 


J*10 4 . 


PS 


(m) 


( 11 ) 


The height above the bed in meters of the CTD (zctd) was then found using the variables 
calculated above. 


z CTD = l-2m + z paro -z r] (m) 


( 12 ) 


2. ADCP 

The ADCP was adjusted to the HAB reference by first determining the height 
offset of the ADCP instrument (1.2m), the thickness of the frame (0.03m), and the 
distance to the first bin (1.61m). By adding each of these components, the HAB of the 
ADCP was approximated to be 3.19m. 

3. BCDV 

Before the HAB for the BCDV could be determined, the range to the time-varying 
bed had to be determined using 30 second backscatter energy. The power of the 
backscatter was computed. The first maximum power above a given threshold was 
selected to be the location of the bed. A threshold of le 5 (in arbitrary units) had to be set 
in order to avoid initial high intensity returns in bubbly conditions. This threshold also 
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extended to the mid ranges of 50cm to 70cm from the instrument where values greater 
than le 5 were sometimes encountered due to a large amount of backscatter from 
suspended sediment. When the threshold was exceeded in this range, the maximum of 
each power computed was determined to be the range to the bed. Ranges to the bottom 
bed bin varied anywhere from 79cm to 81cm over the 2 month observation period. This 
range was then inverted to obtain the HAB for the backscatter and velocity profile 
timeseries. 
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IV. RESULTS 


A. SOLITON OBSERVATIONS 

Internal tidal bores (ITBs) and solitons were observed regularly at the MISO site 
during the two-month experiment. The tidal bores propagated both in a cross-shore (CS) 
and alongshore (AS) direction and contained either single soliton events or packets of 
solitons on the leading edge. These events had a significant impact on the water column 
structure as they passed through the observation site. In the following sections, several 
case studies were chosen to demonstrate and compare the different properties of the CS 
and AS bITBs and their effect on the water column and the bottom boundary layer. 

1. YD 228-228.15 

The internal tide of YD 228-228.15 was a clear example of a soliton that formed 
on the steepening edge of the ITB. Temperature profile timeseries in Figure 14a show 
strong stratification present at the MISO site with the thermocline at 8-10m HAB before 
the arrival of the soliton. 4m isotherm displacements down to 6m HAB occur, causing 
the displaced water to have temperatures 4°C warmer than the background at depth. This 
large displacement and the opposed top layer and bottom layer current in Figure 14b and 
14c indicate that both the soliton and tidal bore are strong mode 1 baroclinic processes. 
The single soliton on the steep edge of the internal bore arrives at the observation site at 
YD 228.08 as seen in Figure 14a and 14b. The duration from start to finish of the soliton 
pulse is approximately lOmin and can be considered as a sharp temporal version of the 
internal tidal bore which has a duration from the leading edge to maximum amplitude of 
approximately 75min. 

Both the ITB and solitons of YD 228 are highly nonlinear baroclinic processes in 
which isotherm displacements and velocities change suddenly with depth and time as 
pulses, rather than sinusoidal oscillations. Temperature profile timeseries show a strong 
near surface stratification present prior to the arrival of the soliton and ITB. As the ITB 
pass in these very shallow conditions, the stratification changes rapidly from near surface 
to near bed. Current vector profiles every lm from 3m to 12m HAB also show a very 
rapid change in currents as they switch sign during the passage of the ITB and soliton in 
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Figure 14. Yearday 228-228.15. (a) Temperature profile timeseries. (b) ADCP CS velocity 
profile timeseries. (c) Current vectors. Vertical alignment of isotherm and current 
displacement are indication of a soliton and internal tidal bore. 
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Figure 14c. Both the ITB and soliton have velocity signatures of 0.10ms" 1 predominantly 
onshore near the surface, indicating that the bore is propagating CS from the shelf break 
directly offshore of the MISO site. Due to the baroclinicity of the approximately two- 
layer system, currents below the bore and soliton are moving with comparable velocities 
of 0.10ms" 1 but in the opposite direction (onshore). The current vectors illustrate how the 
water column is weakly sheared prior to the onset of the bore, with steady currents of 
0.02-0.05ms" 1 moving in a +AS direction near the surface and the bed and with no sudden 
changes occurring from YD 228-228.08. As the soliton approaches the observation site, 
the currents become increasingly offshore from 12m to 8m HAB at YD 228.08, but 
continue to stay rather AS from 7m down to the bed. However, during the passage of the 
soliton and ITB from YD 228.08-228.15, currents at 10m HAB rotate to an onshore 
direction. Within the ITB and soliton, current velocities rotate 180 0 indicating water 
being pushed onshore (-x) direction, while below the soltion and the leading edge of the 
ITB waters are being pulled offshore with the same magnitude. These sudden, pulse-like 
currents associated with solitary internal waves and internal tidal bores have significant 
effects on the water column as they are the dominant currents occurring at the site. 

2. YD 229-229.1 

The packet of solitons occurring during YD 229-229.1 is an example of multiple 
solitons that have formed on the leading edge of an ITB. Temperature profile timeseries 
in Figure 15a shows strong stratification present similar in structure to YD 228, with the 
thermocline at 8-10m HAB prior to the arrival of the bore. Again, approximate 4m 
isotherm displacements occur. These large isotherm displacements and opposed top layer 
and bottom layer current in Figure 15b and c, again show the strong mode 1 baroclinic 
processes controlling these currents. 

Unlike the soliton and ITB found in YD 228-228.15, these 3 packet solitons and 
their ITB are not as steep, but rather weaker and broader, and have decreasing, rank 
ordered amplitude. The duration of these packet solitons are on the order of 15min, 
approximately 1/3 longer than the sharp pulse soliton in YD 228. As in YD 228, 
isotherm displacements and velocities change suddenly with time and depth and occur as 
pulses rather than sinusoidal oscillations. Within the temperature profile timeseries of 
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Figure 15. YD 229-229.1 (a) Temperature profile timeseries. (b) CS ADCP current velocity 
profile timeseries. (c) Current vector timeseries. CS ADCP current velocities indicate 
solitons as part of a wave train. 
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Figure 15 a, strong, near surface stratification is present prior to the arrival of the ITB and 
packet solitons and changes rapidly to near bed as they pass in the shallow conditions. 
Current vector profiles every lm from 3m to 12m HAB show a weaker change down 
through the water column as the packet solitons move through from YD 229.05-229.1. 
However, the CS signatures have approximately the same magnitude as in YD 228 with 
predominately onshore currents of 0.10ms" 1 near the surface and offshore currents of 
comparable velocities below the solitons. As in YD 228, the current vectors in Figure 
15c for YD 229 illustrate a stable water column prior to the arrival of the solitons with 
currents of 0.05ms" 1 in a +AS direction from YD 229-229.05. However, during the 
arrival and passage of the first soliton at YD 229.05, currents from 6-10m HAB rotate to 
an onshore direction. This onshore current decreases in magnitude and becomes more 
+AS as each successive soliton arrives at the observation site. The current vectors 
indicate that within each soliton, water is being pushed onshore, becoming more +AS 
through the packet, while below the soliton, waters are continuously being pulled 
offshore with equal or greater magnitudes. The soliton packet of YD 229 are sudden 
pulse-like currents, but of longer duration and weaker magnitude compared to YD 228 
events, with similar effects on the water column. 

B. INTERNAL TIDE BORE (ITB) OBSERVATIONS 

By expanding the timeseries to a 24-hour period, the effects of the semidiurnal 
ITB over two cycles can be seen more clearly in YDS 228, 229, and 232 in Figures 16, 
17, and 18, respectively. These three one day timeseries illustrate ITBs with decreasing 
steepness, and consequently fewer solitary internal waves (SIWs), over a four day 
interval. 

1. YD 228 

The steep edge of the ITB is clearly seen in YD 228 in the temperature and CS 
velocity profile timeseries of Figure 16a and 16b. Because of missing temperature data, 
the semidiurnal nature of the bore is not seen in the temperature profile timeseries, but the 
stronger first bore can be seen in the CS velocity profile. The first ITB that occurs at YD 
228.1 is much steeper and stronger than the second ITB occurring later on that day as 
seen in the temperature and CS profile timeseries of Figure 16a and 16b. The steep 

leading edge of this first dominant bore enables nonlinearites to develop and form the 
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sharp soliton pulse seen in the previous case study of YD 228-228.15 in Figure 14. The 
steep edge duration of this ITB is on the order of 75min, while the weaker second bore is 
difficult to define. However, both of these tidal bores have a nominal periodicity of 12 
hours. These top and bottom layer opposed currents are not as dominant for the weaker 
and less steep second bore. 

A fourth order low pass Butterworth filter with a cutoff of 5min was used to filter 
out higher frequency solitary internal wave events in order to focus on the effects of the 
ITB on currents throughout the depth of the water column in Figure 16c. Low passed 
filtered current vectors indicate that just prior to the arrival of the ITB from YD 228- 
228.08, currents of O.Olms^-O.OSms" 1 are present in the +AS direction. However, upon 
the arrival of the ITB at YD 228.08, near bed currents rotate onshore with speeds of 0.02- 
0.05ms 1 within the steep leading edge and offshore with equal magnitude just below it. 
As the tidal bore passes the observation site, baroclinic processes push water onshore 
from the surface down to 6m HAB and offshore from 6m HAB to the bottom during the 
duration of the bore from YD 228.08-228.15. At the maximum amplitude displacement 
of the ITB (YD 228.15), currents return to a depth dependent + or - AS direction. 
Because the second bore arriving at YD 228.6 is not as steep and much weaker than the 
first bore, no onshore movement of water is observed within the leading edge, however, a 
slight offshore current is seen directly below it from 5m3m HAB at YD 228.7. +AS 
currents are more prevalent in this weaker bore, possibly indicating that it was not 
generated directly offshore at the shelf break, but rather at the Monterey Canyon and has 
propagated down the bay. 

2. YD 229 

YD 229 is another example very similar to that of YD 228. Again, two ITBs are 
readily observed in the temperature profile timeseries of Figure 17a, with the first one 
being much steeper, with displacements of 8m and currents >0.1 ms" 1 . The first ITB 
arrived at the observation site at YD 229.08 and contains the packet of solitons that was 
discussed in case study 229-229.1. Like the solitons, this bore is not as steep as the first 
bore that occurred during YD 228. The steep edge duration of this first bore is 
approximately 3.5 hours long, which is more than twice as long as the previous example. 
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Figure 16. YD 228-229. (a) Temperature profile timeseries. (b) CS ADCP velocity profile 
timeseries (c) Low pass filtered current vectors. Although the temperature timeseries has 
missing data, the semidiurnal internal tidal bore is evident in the CS ADCP velocity and 
current vector timeseries. 
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The second bore occurs at YD 229.75 and is not as clearly defined in the temperature 
profile timeseries or the velocity profile timeseries, making it difficult to approximate its 
steep edge duration. 

By examining both the temperature and velocity profile timeseries (Figures 17a 
and 17b), the first ITB has large isotherm displacements down to 2m HAB and rapid 
velocity changes with depth and time. Stratification again changes from near surface to 
near bed during the passage of this ITB and CS velocities are observed to be onshore 
within the leading edge and offshore just below it. In contrast, the second ITB has very 
small isotherm displacements of 4m, and no significant changes in the CS velocities. 
Stratification is seen to change only slightly and remains predominantly near surface, 
while no sudden changes in CS velocities are observed between the leading edge and the 
currents below it. 

Figure 17c shows the low passed filtered current vectors for YD 229. Prior to the 
arrival of the first ITB and its soliton packet, current are 0.05 ms 1 in the +AS direction. 
As the soliton packet arrives at YD 229.08, currents from 12m down to 7m HAB rotate 
slightly to an onshore direction, but with a strong +AS component, while the currents 
directly below the solitons, nearer the bed from 5m3m HAB, rotate to an offshore 
direction. In contrast, currents near the peak displacement of the ITB rotate back to a 
+AS direction, but remain offshore beneath it. After the maximum amplitude 
displacement of the ITB, currents reverse to a -AS direction throughout the depth of the 
water column. The strong AS currents that occur in conjunction with the CS currents 
indicate that the ITB has strong vortical mode. The second ITB, in contrast, has no 
significant impact on the currents within the water column as they remain fairly +AS 
throughout its passage. However, a small offshore current is observed under the broad, 
weak leading edge of the bore from 7m-3m HAB at YD 228.85. Once the maximum 
amplitude displacement occurs at YD 229.9, currents return to a +AS direction, in 
contrast to the -AS direction of the first ITB. 
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Figure 17. YD 229-230(a) Temperature profile timeseries. (b) CS ADCP velocity profile 
timeseries (c) Low pass filtered current vectors. Although the temperature timeseries has 
missing data, the semidiurnal internal tidal bore is evident in the CS ADCP velocity and 
vector timeseries. 
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3. 


YD 232 


YD 232 had a contrasting example of two low slope ITBs (Figures 18a and 18b). 
The first bore occurring at YD 232.2 has decreased in steepness when compared to the 
ITBs of YD 228 and 229. The broad, leading edge of this bore is not capable of 
generating any distinct or dominant solitons as in the previous two cases. The second 
bore observed at YD 232.75, like the first, has a broad leading edge that has not 
generated any solitons. However, both of the ITBs occurred with equal magnitude, 
unlike the previous cases where the first bore was significantly dominant over the second 
one. The semidiurnal nature of these two bores is in the CS profile velocity timeseries 
with a period of 12 hours. The steep edge duration of these two ITBs is approximately 4 
hours long, more than twice as long as the steeper edge bores of YDS 228 and 229. 

Mode 1 baroclinic processes are again evident for this case with isotherm 
displacements down to 3m HAB occuring on the leading edge of the two ITBs. During 
the passage of each bore, stratification changes from near surface to near bed as was 
observed in the previous two examples. Changes in CS velocity between the upper and 
lower layers are observed as well, but are much weaker when compared to the steeper 
edge bores of YDS 228 and 229. 

Currents prior to the arrival of each bore are in a +AS direction (Figure 18c). A 
strong vortical mode is present in the currents upon the arrival of the first ITB as currents 
are transitioning from +AS to -AS seen from YD 232.0-232.3. Within this time period, 
currents from 12m5m HAB rotate rapidly from +AS to onshore to -AS. These onshore 
currents are occurring within the leading edge of the ITB. However, just below the 
leading edge of the bore from 5m-3m HAB, currents rotate rapidly from +AS to offshore 
to -AS. Like YD 229, the second bore remains fairly +AS through its passage from YD 
232.6-232.9 with a slight offshore current directly beneath the leading edge from 7m-3m 
HAB occurring at YD 232.8, indicating a weaker, but present vortical mode. Because 
this bore is propagating in a +AS direction, it was most likely generated at the axis of the 
Monterey Canyon to the north of the observation site and was refracted down the bay. 


42 



Temperature Timeseries 


.UvA^M, 



Current Vectors 



0.2 ms 


Figure 18. YD 232-233. (a) Temperature profile timeseries (b) CS ADCP velocity profile 
timeseries shows clear indication of the semidiurnal internal tidal bore, (c) Low pass 
filtered CS vectors illustrate the baroclinicity of the water column. 
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The second bores of YDS 228, 229 and 232 are clear examples of bores that were 
not generated at the shelf break directly offshore, but more likely at the Monterey Canyon 
then subsequently refracted to propagate AS down the bay. Because these bores traveled 
a greater distance to reach the observation site is a possible indication that energy was 
lost due to dissipation as was indicated by their broader leading edges and weaker CS 
velocity signatures. Unlike the first ITBs of YDS 228 and 229, there were no dominant 
soliton events on the leading edge of these secondary bores. The energy loss most likely 
prevented solitons from evolving or caused them to decay before reaching the 
observation site. 

The above case studies indicate that tidal bores moved regularly in both the CS 
and AS direction. The bores generated at the shelf break directly offshore of MISO 
propagated in a CS direction and were observed to contain single or multiple soliton 
depending on the steepness of the leading edge. These events were capable of pushing 
currents onshore near the surface and offshore near the bed. However, the bores 
generated at the axis of the Monterey Canyon to the north of MISO, propagated in an AS 
direction and were much broader and weaker and therefore, contained no solitons on the 
leading edge. 

4. Tidal Cycle 

The ITBs observed over the four-day period from YDS 228-232 are not closely 
tied to the barotropic tidal cycle. The ITBs pass through the observation site at different 
times each day. The first ITB of YD 228 occurs near the spring tide of the semidiurnal 
barotropic tide (Figure 19). The bore happens to pass through the observation site during 
the ebb cycle of the first high tide, at a time determined by the propagation delay between 
the generation and observation sites (Figure 20). This bore was generated during the 
transition from ebb to flood flow of the previous 2 nd high tide (YD 227). However, the 
second, weaker bore passes through the observation site during the ebb cycle of the 
second high tide and was generated at the shelf break during the transition from ebb to 
flow of the previous first high tide (YD 228). This same pattern is observed for the ITBs 
of YDS 229 and 232 in Figure 20, where the steeper first bore (large star) passes the 
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Figure 19. Tidal timeseries from YD 184-244. Periods of maximum amplitude are known as 
spring tides and period of minimum amplitude are neap tides. YD 228 is occurring near 
the spring tide of the tidal cycle. 


observation site during the ebb cycle of the first high tide, being generated within the 
second high tide of the previous day, while the second weaker bore (small star) passes 
during the ebb cycle of the second high tide, being generated within the first high tide of 
that day. The amplitude displacement that occurs between max flood and max ebb for 
the second high tide is 5m which is much smaller than the 20m displacement for the first 
high tide. Less energy is lost for the smaller displacement allowing the bores generated 
from it to be the stonger and steeper. However, because more energy is lost for larger 
displacements, the ITBs generated from it tend to be weaker and broader. The ITBs are 
also assisted by the respective tide that carries them through the observation site. The 
strong ITBs pass through MISO during the larger magnitude first high tide, while the 
weak ITBs pass through during the lower magnitude second high tide. 
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Figure 20. Tidal Pressure Timeseries showing the four day period between YDS 228-232. 
Vertical axis is in m (HAB) x 10. ITBs are indicated in the figure as stars. The steeper 
first ITBs (large stars) that were observed at MISO occurred during the ebb flow of the 
first high tide, while the weaker second ITBs (small stars) occurred during the ebb flow 
of the second high tide. 
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C. BOTTOM BOUNDARY LAYER OBSERVATIONS 

The internal tidal bores and accompanying solitons at the MISO site produced 
currents near the bed significantly greater than background levels. The pulses of current 
generated a deep boundary layer above the bed, and it is of interest to see whether these 
currents contributed significantly to sediment suspension. The evolution of a logarithmic 
layer was observed as the ITB and soliton currents moved over the sandy, mobile bed. 
This resulted in a wall bounded shear flow. By using the “Law of the Wall'’ in the region 
far from the viscous sublayer, a dimensional analysis yields (Kundu 1990) 


dU _ u* 
dz kz 

Integration gives 


(13) 


U(z) = — fo(z) +const. (14) 

k 

A profile plot using log(z) versus current speed, reveals a logarithmic layer as a straight 
line. In this region, U (current speed) is related to u* by the equation 


and so 


U(z) = -H-) 


(15) 


U(z)k 
ln( z) 


(16) 


where 


U(z) 
ln( z) 


is the inverse slope of the logarithmic layer, and k is the Von Karman 


constant, approximately 0.41. Values of u* under the SIWs and ITBs were on the order 
of 0.01-0.03ms" 1 as ITBs passed the MISO site. These high friction velocity values 
indicate that significant stress occurs due to friction drag as the SIW flow field interacts 
with the dynamically rough bed. Because u* values scale roughly with the mean current, 
as the mean current approaches zero, so does u*, and the friction velocity becomes 
difficult to estimate from fits to the log profiles. Once the friction velocity has been 
estimated, the shear stress above the bed is then estimated from 
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(17) 


T = p M* 

where p is the density of the water. 

Very near the bed layer (Omm) in the viscous sub-layer, molecular viscosity 
effects dominate fluid movement. This region is too small a scale to be studied in our 
bottom boundary analysis. In coastal areas with sandy beds, ripples have been observed 
to dominate the apparent roughness over mobile, sandy beds. Since surface gravity 
waves result in oscillatory fluid motion down to the bed, a thin wave boundary layer 
forms, which tends to decouple fluid movement from above. This occurs when there are 
significant levels of long period surface gravity wave energy within the water column. 
This wave boundary layer causes the mean currents to attenuate very slowly above the 
bed in the log layer (Figure 21). Without the presence of this wave boundary layer, the 
logarithmic layer has a direct interception to zero current (Figure 22). 

A logarithmic layer evolves above the bed as the baroclinic flows from SIWs and 
ITBs pass the MISO site. This turbulent shear layer is bounded above by strong 
stratification associated with baroclinic events. Because this analysis focuses on the 
effects of solitons on the bottom boundary layer, stress timeseries are estimated from the 
near-bed logarithmic layer. Within the logarithmic layer, there is an inertial transfer from 
the viscous sublayer to the wave boundary layer and then the logarithmic outer layer by 
inviscid nonlinear processes. 

Polynomial fitting and regression techniques have been used to automatically 
determine the slope of this logarithmic layer to produce timeseries of u*. Between 4cm 
above the bed to 7cm, a linear regression line was calculated and its regression 
coefficient was found. Points were added one cm at a time until the correlation began to 
decrease. The fit that scored the highest correlation value above 0.80 was used to 
determine the slope of the logarithmic layer (Figure 23). This slope was then used to find 
the zero velocity intercept, giving a coarse estimate of the apparent roughness of the bed. 
These values of Z 0 were only resolved to 1cm, but provided Z 0 estimates for high 
roughness conditions. With steady flow over the bed, a wake is generated behind each 
bump in the bed and stress is transmitted to the bed by way of a pressure drag on the 
roughness elements. 


48 



10 


Regression of Log Layer r -0 9275 Day 206.1214 


*3ep 
* bc^y 


10 


£ 


10 


0 


f 




10 


-01 


-0.05 


005 


Vector Vetac *y m/s 


—L_ 

01 


Figure 21. A plot of total velocity with logio(z) vertical coordinate for a 30sec mean current. 
Slow attenuation of the logarithmic layer for YD 206.12, during a high energy wave day 
caused by the wave boundary layer. 
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Figure 22. A plot of total velocity with logio(z) vertical coordinate for a 30sec mean current. 
Fast attenuation of the logarithmic layer for YD 219.2, during a low energy wave day 
caused by the absence of the wave boundary layer. 
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Figure 23. Linear regression of the logarithmic layer for YD 228.08 and calculation of the 
slope. 


The effects of SIWs, ITBs, and surface gravity waves on the bottom boundary 
layer and sediment suspension are considered over a range of forcing conditions in 
Figures 24-29. In each of these 3-hour timeseries, broken into two groups of high and 
low surface gravity wave energy days, the RMS wave velocity at 3m HAB, 3m HAB 
mean ADCP and 4-18cm HAB mean BCDV current components, Shields Parameter 
(Chapter V), relative backscatter profiles (that scale with sediment concentration), and 
friction velocity are shown. 

1. High Wave Energy 

Relatively high wave energy days are shown in Figures 24, 25 and 26 for YDS 
200.8, 206.05, and 219.05. YD 200.8 shown in Figure 24 was one of the first high wave 
energy days encountered during the two month observation period. During this 3-hour 

period RMS wave velocities were approximately 0.12ms" 1 as seen in Figure 24a. Effects 
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of this high wave energy, due primarily to surface gravity wave swell, is seen in the 
backscatter profile where shades of white indicate high backscatter and shades of gray 
and black indicate areas of low backscatter (Figure 24d). These backscatter values scale 
with sediment concentration. The increase in backscatter from black to white values 
corresponds to sediment concentrations that have increased by 3 orders of magnitude. 
Figure 24b shows the 3m HAB mean ADCP (bold line) and 4-18cm mean BCDV CS 
velocities. A sharp soliton pulse in the lower 40cm with an onshore velocity of 0.10ms" 1 
is observed at YD 200.92 on the leading edge of an ITB with a CS velocity of the same 
magnitude. The presence of the strong swell wave boundary layer has also helped to 
decouple the fluid from above resulting in a small vertical velocity gradient from 3m 
HAB down to the lower 40cm. 

The instantaneous Shields Parameter (O'), which is an indication of when 
sediment motion is initiated (to be discussed further in Chapter V), is shown in Figure 
24c. After a sensitivity analysis, twice Jonsonn’s (1967) wave friction factor was used to 
illustrate the effect of high surface gravity wave forcing. The critical Shields Number 
(0 cnr ) for the well-sorted sediment at the site is 0.05, and is a rough indicator of when 
sediment suspension will occur. Prior to the arrival of the soliton and ITB at the 
observation site, 0' has a value of 0.06-0.08, indicating that the surface gravity wave 
forcing is strong enough to suspend sediment past 0 cnt . However, during the arrival and 
passage of the soliton and ITB, 0' increases significantly to values of 0.06-0.12, well 
above 0 rnt and the values that were generated by the surface gravity waves, indicating 

that the soliton and ITB coupled with surface gravity wave forcing provided more than 
enough stress to suspend sediment. It is evident, therefore, that the stress generated by 
high surface gravity wave forcing results in wave-dominated friction velocity values of 
1.5-2.0cms" 1 from YD 200.8-200.9 in Figure 24e. The friction velocity is seen to 
increase significantly when a SIW and ITB arrive with strong near bed currents. The 
bold line indicates the critical friction velocity of 1.5 lcms" 1 (to be discussed further in 
Chapter V), in which values exceeding this will cause sediment suspension. The high 
suspended sediment profiles (Figure 24d) show the surface gravity wave forcing to be 
sufficient to produce a large suspended load. 
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Figure 24. YD 200.8-201. (a) RMS Wave Velocity (b) CS ADCP (3m HAB) (bold line) and 
BCDV (4-18cm HAB) (c) Instantaneous Shields Parameter 0 mf =O.O5 is denoted by the 

threshold line (d) Backscatter profile timeseries (e) Friction Velocity u* u* cnt =\ .5 lcms -1 is 
denoted by the threshold line. 
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YD 206.05-206.2 in Figure 25 is a moderate surface gravity wave energy day 
with RMS wave velocities of 0.07 ms" 1 at 3m HAB (Figure 25 a). The backscatter profile 
timeseries of Figure 25d shows an increase in sediment suspension that occurs when an 
ITB with multiple solitons passes through the observation site at YD 206.12. Unlike the 
previous example, where high backscatter was occurring from strong surface gravity 
wave forcing, the SIW and ITBs contribute significantly to suspension here. In this case, 
the surface gravity waves are not as strong in causing sediment suspension. Rather it is 
the ITB and solitons coupled with the moderate surface gravity wave energy that is able 
to generate suspended sediment. Figure 25b reveals the multiple weak soliton packets, 
which are both weaker and less steep, than the previous example shown. The peaks of 
the weak soliton currents at 3m HAB are nearly out of phase near the bed. The effects of 
these weaker events are evident in the 0'and u* timeseries of Figure 25c and 25e as they 
only help to slightly increase the values past the critical thresholds. However, the 
backscatter timeseries does indicate that these events, though weak, when coupled with a 
moderate surface gravity wave forcing, are capable of suspending sediment in the water 
column. 

YD 219.05 in Figure 26 is another example of a high surface gravity wave energy 
day, but coupled with a rather weak soliton and ITB. RMS wave velocities shown in 
Figure 26a are high, exceeding 0.12ms" 1 . This high surface gravity wave forcing is again 
evident in the backscatter profile timeseries of Figure 26d, where the lighter areas are 
approximately 2-3 orders greater in sediment concentration than the darker areas. The 
highest backscatter occurs at the peak RMS wave velocity of 0.12ms" 1 at YD 219.07, 
indicating that surface gravity waves are causing sediment to become suspended above 
the bed. As the RMS wave velocity slowly decreases within the 3-hour timeseries to 
0.08ms" 1 , high backscatter still occurs. This is assisted by the arrival of solitons and an 
ITB at YD 219.12. Multiple solitons are observed in the CS ADCP and BCDV velocities 
of Figure 26b with rather weak near bed velocities of 0.02-0.06ms" 1 . As in YD 200.8, the 
wave boundary layer reduces vertical shear in CS velocities between 3cm and 40cm. 
However, unlike in YD 200.8, these solitons and the ITB appear broader and not as steep 
and therefore, do not have as significant of an impact on increasing values ofO'and u* as 
seen in Figure 26c and 26e. These two timeseries do indicate that the strong surface 
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Figure 25. YD 206.05-206.2. (a) RMS Wave Velocity (b) CS ADCP(3m HAB) (bold line) 
and BCDV (4-18cm HAB) (c) Instantaneous Shields Parameter 0 crif =O.O5 is denoted by 
the threshold line (d) Backscatter profile timeseries (e) Friction Velocity u* 
u*criF^ 1 • 51 eras 1 is denoted by the threshold line. 
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gravity waves are suspending sediment, but that the solitons and ITB help to sustain 
sediment suspension as the wave energy decreases. 

2. Low Wave Energy 

YDS 191, 214, and 228 (Figures 27, 28 and 29) are counter examples to those 
described above. YD 191 is a clear case of low surface gravity wave energy coupled 
with strong soliton pulses and a steep edge ITB. Here, the low surface gravity wave 
energy in the water column is reflected in the RMS current velocities of 0.04 ms" 1 (Figure 
27a). Because the wave boundary layer is much thinner and weaker at these low values, 
CS ADCP and BCDV velocity profiles of Figure 27b indicate a significant vertical 
velocity gradient in the lower 40cm. The near bed velocities also indicate the presence of 
strong solitons and a steep ITB with onshore velocities of 0.10-0.12ms" 1 , which result in 
large suspended sediment loads seen in the backscatter profile of Figure 27d. At the 
maximum amplitude displacement of the soliton during YD 191.16, a significant increase 
in backscatter is seen as well as during the multiple soliton events on the ITB occurring 
during YD 191.24-191.26. Unlike during the high wave energy days, where it was 
shown that surface gravity wave forcing was strong enough to suspend sediment, it is the 
solitons and the ITB that are generating the friction velocities seen in Figure 27d that lead 
to sediment suspension. The sharp pulses of the solitons in this case, are seen to only 
slightly increase values of 0' and u* past the critical values, but are still capable of 
causing sediment suspension as indicated in the backscatter timeseries. 

YD 214.27 in Figure 28 is an example of a low wave energy day coupled with a 
very broad ITB at the beginning of the 3-hour time period and distinct soltion pulses and 
a second ITB towards the end of the period. The low wave energy is again evident in the 
RMS wave velocities of 0.04ms" 1 in Figure 28a. Despite the low surface gravity wave 
forcing, the backscatter timeseries indicates that sediment suspension is occurring (Figure 
28d). When compared with the CS ADCP and BCDV velocity timeseries, the increases 
in backscatter occur during the passage of the SIWs and ITBs. The near bed CS 
velocities of Figure 28b indicate multiple SIWs occuring at YD 214.22 and YD214.34 
with strong offshore velocities of 0.07-0.10ms" 1 . These ITBs are strong enough to 
generate friction velocities capable of suspending sediment as indicated by the 0' and u* 
timeseries of Figure 28c and 28d. 
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Figure 26. YD 219.05-219.2. (a) RMS Wave Velocity (b) CS ADCP (3m HAB) (bold line) 
and BCDV (4-18cm HAB) (c) Instantaneous Shields Parameter 0 crif =O.O5 is denoted by 

the threshold line (d) Backscatter profile timeseries (e) Friction Velocity u* 
u* crit =\ .5 lcms" 1 is denoted by the threshold line. 
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Figure 27. YD 191.1-191.3. (a) RMS Wave Velocity (b) CS ADCP (3m HAB) (bold line) 
and BCDV (4-18cm HAB) (c) Instantaneous Shields Parameter Q crit =0.05 is denoted by 

the threshold line (d) Backscatter profile timeseries (e) Friction Velocity u* 
u* cnt =\ .5 lcms -1 is denoted by the threshold line. 


57 




















RMS Wave Velocity 



</) 

E 0.02 


CS ADCP (3m HAB) BCDV Velocity 



214,2 214.22 214 24 214 26 214 28 214 3 214 32 214 34 214 36 214 38 214 4 

yearday 2002 


Figure 28. YD 214.2-214.4. (a) RMS Wave Velocity (b) CS ADCP (3m HAB) (bold line) 
and BCDV (4-18cm HAB) (c) Instantaneous Shields Parameter 0 crif =O.O5 is denoted by 
the threshold line (d) Backscatter profile timeseries (e) Friction Velocity u* 
u*criF^ 1 • 51 eras 1 is denoted by the threshold line. 
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YD 228-228.15 is another low surface gravity wave energy coupled with a single 
soliton and a steep edged ITB shown in Figure 29. The low wave energy is again evident 
in the RMS wave velocities of Figure 29a where values are on the order of 0.4ms" 1 . This 
low surface gravity wave forcing is clearly reflected in the backscatter timeseries of 
Figure 29d, where shades of black and gray completely dominate the timeseries, 
indicating no sediment suspension is occurring above the bed. However, light areas 1.5-2 
orders greater in magnitude are occurring at approximately YD 228.08. By comparing 
this to the near bed CS ADCP and BCDV velocities in Figure 28b, it is seen that the 
single soliton and steep leading edge of the ITB are again the main causes in backscatter 
increase. These events are generating weaker near bed currents up to 0.04ms" 1 . 
Timeseries of O' and u * indicate periods that exceed the critical values, although some of 
these values appear to be bad data points. 
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Figure 29. YD 228-228.15. (a) RMS Wave Velocity (b) CS ADCP (3m HAB) (bold line) 
and BCDV (4-18cm HAB) (c) Instantaneous Shields Parameter 0 crit = 0.05 is denoted by 
the threshold line (d) Backscatter profile timeseries (e) Friction Velocity u* 
u* cnt =\ .51cms _l is denoted by the threshold line. This is a low wave energy day, which is 
evident in the backscatter (areas of black/dark gray). Solitons are evident in both the CS 
BCDV and ADCP and are strong enough to cause sediment suspension (areas of white). 
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V. DISCUSSION 


Solitary internal waves (SIWs) and internal tidal bores (ITBs) acting with surface 
gravity wave forcing frequently provide the necessary stress above the bed to suspend 
sediment and organisms that live on the bed into the water column. We hypothesize that 
the strong, unidirectional flows associated with SIWs and ITBs result in significant net 
cross-shelf transport of this suspended load in comparison with background currents at 
this site. This cross-shelf transport of organisms can affect migration, reproduction, and 
feeding of these and other organisms. 

Solitons and ITBs were observed daily during July and August 2002. At least 2 
ITBs occurred daily, with one usually being more dominant than the other and moving in 
either a cross-shore (CS) or alongshore (AS) direction. Multiple solitons were frequently 
observed on the strong CS bores in contrast to the AS bores that carried with them at 
most only one dominant soliton. Because these events occurred during most tidal cycles 
during our observation period, they are potentially responsible for a large part of 
sediment and biologic entrainment and transport within on the inner shelf within 
Monterey Bay. 

These events consist of strong unidirectional pulses of current, unlike the 
oscillatory motion of the waves. Therefore, wave forcing will help to suspend and 
entrain sediment, but solitons and tidal bores provide the dominant forcing for net 
transport. 

A. SEDIMENT SUSPENSION 

In the results section, it was shown that solitons and internal tidal bores interacted 
with the bottom boundary layer through the formation of a logarithmic layer due to skin 
friction drag along the bed and that the flow over the bed caused by these events can be 
considered as a wall bounded shear flow. From the slope of this logarithmic layer, the 
friction velocity was determined and the current stress (x current ) above the bed can be 
estimated as 

^current ~ P M * (18) 
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where p is the density of the water and w* is the friction velocity. 

By combining the stress due to the current and the stress due to the surface gravity 
waves (x wave ), the total bed stress ( x hed ) can be determined. 


T — T -LT 

bed current wave 


(19) 


The stresses from the current and the wave can be broken down into individual 
components for the current and the wave 


^ bed ^ SIW IT LF wave 


( 20 ) 


where x SIW is the soliton component, x /7 is the internal tidal bore component, and x LF is 
the low frequency component consisting of wind and tidal forcing. The SIW and IT 
contributions have been isolated temporarily in this analysis, while x LF is estimated from 
conditions when ITBs and SIWs are absent. 

As stated above, the stress from the current can be directly found through the 
friction velocity that was derived from the slope of the logarithmic layer by using the 
Law of the Wall. The stress from the wave can be further broken down to 
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( 21 ) 


where f w is the wave friction factor following Jonsson (1967). This is obtained using 
the median grain size of the sediment (dso) and the RMS wave orbital diameter (a) that 
can be obtained by integrating the RMS wave velocities (u 2 waveRMS + v 2 vaveRMS ), measured 

over a bandwidth fom 30Hz to the Nyquist frequency (0.2Hz), with respect to their 
frequency, giving 

/ w = exp[5.213(^<L) 0194 - 5.977] (22) 

a 

Inspection of timeseries in low current, high wave conditions suggests a factor of 
2 f w better matches critical Shields number criteria, and this has been used in Figures 24- 
29. 
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Figure 30. The median grain size (dso) at the MISO site was found to be 0.19mm using in- 
situ sediment samples taken from April and June 2002. These samples represent a typical 
log-normal distribution of a single population of well sorted sand in Monterey Bay. 


A median grain size (d^o) of 0.19mm was measured using in-situ sediment 
samples taken in April and June 2002 from the MISO site (Figure 30). 

Once the total stress felt above the bed was computed, it can be determined 
whether or not the stress was sufficient enough to lift and suspend sediment. The point at 
which sediment suspension is initiated occurs when the lifting force due to the stress is 
strong enough to overcome the gravitational restoring forces on the grain. The Shields 
Parameter (0') (Shields 1936) states this simple ratio as 


0 ' = 


T -l-T 

current wave 


PO-1 )gd 


50 


(23) 
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where s is the ratio of sediment density to that of water ( ^‘ ;gdime ” f ) taken to be 2.65 for 

p 

r water 

quartz, g is the gravitational constant at 9.8 ms , and dsc is the median grain size found to 
be 0.19mm. This value represents the instantaneous Shields Number at a given time for 
the given forcing conditions. 

Sediment suspension is initiated when the values for 0' exceed a critical Shields 
number (0 f ) of 0.05, which requires a critical friction velocity u i( of 1.51cms _1 for the 
sediment observed at the MISO site. 

In this analysis, the objective is to determine which factors are responsible for 
sediment suspension near the bed. Solitary internal wave forcing, internal tidal forcing, 
lower frequency forcing such as winds, as well as wave forcing will all be compared to 
determine which contributes the most stress that leads to sediment suspension. 

1. Wind Forcing 

Winds are predominantly from the northwest year round in Monterey Bay except 
when strong storms move through the region and disrupt the pattern. Figure 30 shows the 
hourly sampled winds for Monterey Bay as measured by the MISO MET instrument 
located directly onshore from the MISO site for July and August 2002. These winds have 
been rotated into a CS normal coordinate system and indicate that winds are 
predominantly onshore with a positive AS component. Wind forcing should only be felt 
in the near surface waters since stratification prevents wind driven motion from 
continuing down through the water column and affecting currents near the bed. By 
comparing the CS wind components with the CS current components at 10m HAB and 
3m HAB, the effects of wind forcing on currents near the surface and near the bed can be 
determined. 

A cross correlation between hourly sampled CS winds and CS currents at both 
10m and 3m HAB in Figure 31 shows that there is a higher correlation between wind and 
current at 10m HAB vice 3m HAB for AS components. CS components at 10m HABhad 
a negative correlation while at 3m HAB no significant correlation was sen. The AS 
components at 10m HAB had a positive correlation that was significantly larger than the 
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Figure 31. MISO wind timeseries for YD 184-244. Hourly averaged winds show strong 
diurnal variability and a predominantly northwest duration with a strong onshore 
component. 


correlation at 3m HAB. The low correlations for both AS and CS components indicate 
that the winds within Monterey Bay have only a slight influence on currents 2m below 
the surface and almost no influence on currents 10m below the surface for both CS and 
AS components. This confirms that the stratification present within the water column 
largely decouples wind driven motion the lower part of the water column, and therefore 
prevents wind forcing from having any significant effect on the contribution to the stress 
above the bed. 
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Correlation-Wind and 10m HAB R 2 - 0.18342 



Correlation-Wind and 3m HAB R 2 - 0 012996 



Correlation-Wind and 10m HAB R 2 - 0.092277 



Correlation-Wind and 3m HAB R 2 - 0.0013246 



Figure 32. Correlations between wind and current for YD 200-230. (a) AS wind versus AS 
currents have a higher correlation at 10m HAB vice 3m HAB. (b) CS wind versus CS 
currents follow the same pattern. 
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2. Internal Tidal Bores, Solitons, and Wave Forcing 

The solitons and internal tidal bores observed during the two month observation 
period were frequently the dominant currents felt in the bottom 1/3 of the water column, 
indicating that they most likely disturbed the bed. The resulting current flow over the 
rough bed generates a logarithmic layer typically extending 20-300cm above the bed, 
allowing ITB and SIW friction velocity timeseries to be generated (see Figures 24e-29e). 
With this friction velocity, the corresponding stress above the bed due to the event was 
calculated as shown by equation 17. 

By examining the RMS wave velocities over the same two-month observation 
period, stress produced by waves can be compared to the stress produced by the solitons 
and internal tidal bores. From this, it can be determined whether the dominant stress felt 
above the bed is due to the waves or to the solitons and tidal bores. Figures 24e, 25e, and 
26e show the friction velocity timeseries for the three case studies of high wave energy 
days (YDS 200.8, 206.05, and 219.05). The surface gravity wave velocities as well as 
the solitons and tidal bores that occurred during each day were strong enough to generate 
a logarithmic layer above the bed and, therefore, allow a friction velocity to be estimated. 
RMS velocities due to surface gravity waves ranged from 0.04-0.12 ms 1 causing u* 
values to exceed the critical friction velocity, indicating that strong wave forcing coupled 
with a soliton and internal tidal bore event provided enough friction velocity to suspend 
sediment. Figure 33 summarizes the estimated u* waV e and u* current that comprises the 
SIWs and ITBs. The surface gravity wave energy is shown to have significantly higher 
values of u* when compared to u* values from the baroclinic currents, except for the 
strong ITB events starting on YD 200.92. However, by combining these two u* 
components, the overall u* greatly exceeds the critical value, so the high suspended load 
seen in the acoustic backscatter profiles in Figures 24d-26d are expected. For these 
forcing conditions, surface gravity wave forcing has proved itself to be the dominant 
factor in contributing to bed stress with the solitons and tidal bores enhancing the stress 
above the bed. 
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Figure 33. Friction velocity timeseries ( u *) for high wave energy days (YDS 200.8, 206.05, 
and 219.05). u* due to the waves is significantly higher than u* due to the solitons and 
internal tidal bores. In these cases, surface gravity wave energy is suspending sediment 
inferred from the acoustic backscatter. 


The counter examples of low wave energy days (YDS 191.1, 214.2, and 228) are 
shown in Figures 25d-28d respectively. In these cases, the surface gravity wave energy 
was not strong enough to cause sediment suspension by itself, but rather dependent on the 
additional bed stress due to the solitons and internal tidal bores that occurred during each 
time period. The total friction velocities generated by the waves and currents were below 
the critical value of 1.51cms _1 . Figure 34 shows the friction \elocity due to the waves 
and the currents. Unlike the high wave energy cases where u* waves was significantly 
higher than u* current , the low wave energy days cause u* waV es to be on the same order or 
lower than u* cur rent- 
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Figure 34. Friction velocity timeseries (w*) for low wave energy days (YDS 191.1, 214.2, 
and 228.0). u* due to the waves is significantly lower than u* due to the solitons and 
internal tidal bores. In these cases, surface gravity wave energy is not suspending 
sediment inferred from the acoustic backscatter. 


The different u* values for high and low surface gravity wave forcing indicates 
that periods of low wave forcing play no significant role in contributing to bed stress and 
sediment suspension. The main contribution to bed stress for low surface gravity wave 
forcing was due primarily to solitons and ITBs. However, strong wave forcing creates a 
wave boundary layer in the lower portion of the logarithmic layer that decouples the fluid 
from above and creates more stress above the bed. In a log profile of current velocity 
versus depth, a sharp bend in the log layer is seen in the lower 1-4 cm (Figure 21) as well 
as higher velocities up through the layer. When coupled with a soliton or internal tidal 
bore event, the stress above the bed was seen to increase significantly, resulting in 
sediment suspension. 
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B. NET TRANSPORT 


Unlike solitons and internal tidal bores which are strong unidirectional pulses of 
current, oscillatory motion due to surface gravity waves are largely linear and sinusoidal 
which over time, averages to no overall net transport. We hypothesize that the net 
sediment displacement and water column transport at this site is the dominated by 
currents induced by solitons and ITBs. The net transport of water through the water 
column was analyzed by integrating 10m and 3m HAB ADCP velocities along with 
30cm, 10cm and 3cm HAB BCDV velocities to obtain water column displacement at 
each level for selected time intervals. The six case studies of high and low wave energy 
events will be used again to analyze the net transport associated with each event. 

For the high wave energy days (YDS 200.8, 206.05, and 219.05) in Figures 24-26 
where a strong onshore tidal bore and soltions were observed, net transport in Figure 35 
is shown by integrating over the period of the soliton and internal tidal bore at selected 
HABs. For the YD 200.92 case in Figure 35 a, net transport is directed offshore and 
towards Monterey due to the pulse like current of the soliton and ITB (see timeseries in 
Figure 24b). However, at 10m HAB, the water has been displaced in an onshore 
direction as is characteristic with CS bores and solitons. In the YD 206.05 case (see 
timeseries in Figure 25b), displacement at 10m HAB has a strong AS component in 
combination with its CS movement resulting in a loop-like motion that is likely the result 
of a more vortical mode baroclinic bore passing the observation site (Figure 35b). YD 
219.05 (see timeseries in Figure 26b) shows quite the opposite results in net transport. In 
this case, water from 10m down is being displaced in an offshore direction. However, as 
the bore continues to pass through, water near the bed undergoes several CS oscillations 
associated with more linear near-bed internal waves (Figure 35c). 

The same integration was done for the case studies of low wave energy days. For 
YD 191.1 (see timeseries in Figure 27b), waters below 3m HAB are being transported 
offshore during the entire length of the soliton and bore event. At 3m HAB and 10m 
HAB, there is dominantly onshore transport with a depth dependent AS component 
(Figure36a). For YD 214.2 (see timeseries in Figure 28b), net transport for levels 3m to 
3cm HAB is being directed offshore with a nearly equal -AS component (Figure 36b) 
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and a less clear baroclinic structure, since transport at all the selected levels is largely in 
the same direction. YD 228 (see timeseries in Figure 29b) was described as an example 
of a strong CS propagating bore. This CS movement is seen in the displacement at all 
levels. Onshore transport is occurring at 10m HAB, while from 3m-3cm HAB, a net 
offshore transport is occurring within the water column (Figure 36c). 

It has been shown that significant wave forcing, along with solitons and ITBs can 
generate the required friction velocity to overcome the u* cr it and provide the necessary 
stress to initiate sediment motion past the critical Shields Number. Once sediment has 
been suspended and entrained in different levels above the bed, the strong unidirectional 
pulses of solitons and tidal bores will cause a significant net transport within the water 
column. For ITBs generated at the shelf break directly offshore of MISO, net transport is 
directed largely cross-shelf. However, for those bores generated at the axis of the 
Monterey Canyon to the north, net transport is directed more AS, but with a more 
vortical, looped motion for the less steep ITBs. 

From these observations, we can conclude that surface gravity wave forcing is 
significant in suspending sediment during high wave energy days, but rather insignificant 
for low wave energy days. It was also shown that stronger solitons and tidal bores 
produced sufficient bed stress to suspend sediment with or without wave forcing. Time 
integration of velocity profiles show that SIWs and ITBs were largely responsible for the 
net transport of water and sediment within the water column. 
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Day 206.1-206.2 



Day 210.1 -210 2 



Figure 35. Time integrations for strong tidal bore event and soliton events occurring during 
YDS 200.8 (a), 206.05 (b), and 219.05 (c) showing net transport at different levels within 
the water column. 
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Figure 36. Time integrations for strong tidal bore event and soliton events occurring during 
YDS 191.1 (a), 214.2 (b), and 228.0 (c) showing net transport at different levels within 
the water column. 
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CONCLUSIONS 


Strong internal tidal bores and solitons were observed very close to the nearshore, 
less than 600 km from the surfzone at the inner shelf MISO site. These highly baroclinic 
events were evident in both the temperature and velocity profile timeseries in which 
stratification changed from near surface to near bed and velocities changed from onshore 
at the surface to offshore near the bed. These events extended down to the bottom 1/3 of 
the water column as indicated by the velocity and backscatter profile timeseries. 

Internal tidal bores were observed daily moving in either a CS or AS direction. 
The AS bores were most likely generated at the axis of the Monterey Canyon and were 
refracted to the southern end of the bay. The leading edge of these bores did not 
degenerate into multiple solitons as their CS counterparts did. The CS bores were 
generated offshore at the continental shelf break and were not as frequent as the AS 
bores. However, the leading edges of these bores were able to produce multiple solitons. 

Using wall bounded shear flow theory, it has been demonstrated that a 
logarithmic layer forms due to friction drag as currents flow over the bottom boundary. 
From this logarithmic layer, the friction velocity and associated stresses were determined 
for the soliton and tidal bore events and compared with the stresses produced by wave 
forcing. 

The analysis of sediment suspension and transport under the combined forcing of 
solitary internal waves, internal tidal bores, and surface gravity waves at the MISO site 
has shown that waves are the dominant force to suspend sediment. However, solitons 
and internal tidal bores were found to be the dominant force to transport sediment near 
the bed. 
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